Intra- valley Spin-triplet p + ip Superconducting Pairing in Lightly Doped Graphene 
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We analyze various possible superconducting pairing states and their relative stabilities in lightly 
doped graphene. We show that, when inter-sublattice electron-electron attractive interaction dom- 
inates and Fermi level is close to Dirac points, the system will favor intra-valley spin-triplet p + ip 
pairing state. Based on the novel pairing state, we further propose a scheme for doing topological 
quantum computation in graphene by engineering local strain fields and external magnetic fields. 
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Graphene has attracted great experimental and theo- 
retical interests since its successful fabrication. 0,0 The 
unique electronic structure of graphene, characterized by 
massless relativistic particle-like dispersion near two in- 
equivalent corner points of Brillouin zone, gives rise to 
exotic physical properties. Q One of interesting ques- 
tions would be: how electrons could be paired in such a 
system? There had been many theoretical studies that 
explore the possibility of the superconductivity in sin- 
gle/multilayer graphenes. In particular, Ref. predicts a 
spin singlet p+ip pairing state in the presence of effective 
attractive e-e interaction between the nearest neighbor- 
ing sites. More exotic possibilities such as chiral s-wave, 
d-wave and /-wave pairings are also predicted. First 
principles calculation has suggested that superconductiv- 
ity could be induced by doping alkaline adatoms. Re- 
cently, superconductivity has been realized in potassium- 
doped few-layer graphene. 

In this Brief Report, we investigate possibility of ob- 
taining intra-valley spin-triplet p + ip superconducting 
state (see below) in lightly doped graphene. Unlike the 
spin-singlet p + ip pairing state, Q the intra-valley spin 
triplet p + ip state may support topological excitations of 
Major ana fermions, and could be useful for implement- 
ing topological quantum computation. We carry out a 
systematic investigation on the possible superconductiv- 
ity pairing states and their relative stabilities in lightly 
doped graphene, taking account of multiple degree of 
freedoms of orbit, spin and valley. We find that, when 
effective e — e attractive interaction between the nearest 
neighboring sites (inter-sublattice) dominates and Fermi 
level is close to the Dirac points, the intra-valley spin- 
triplet p + ip pairing state will be favored. Based upon 
the novel pairing state, we propose a scheme for doing 
quantum computation in graphene by engineering local 
strain fields and external magnetic fields. 



We start our investigation from a simplified lattice 
model introduced in Ref. 0, which assumes that effective 
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Figure 1: Two scenarios of pairings: (a) intervalley pairing, 
(b) intravalley pairing. Qo = (0,47r/(3V3a)). 



e — e interaction has the form: 
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where go and g\ denote the strength of onsite interaction 
and interaction between the nearest neighboring (NN) 
sites, respectively. a is (b\ s ) and ai s (bi S ) are the creation 
and annihilation operators on the unit cell i with the spin 
index s =t, 4-j an d in the sublattice A(B). denotes 
all pairs of NN sites. 

We rewrite Eq. |I| in the basis of non-interacting 
eigenstates: rji(k) = (f/\/2) [f, — Zcxp (— 'upk)] T , where 
/ = f (— f) for conduction (valance) band, exp[iip (k)} = 
Tfc/lTfc|, Ik = T, Si eikSi for *i = a(3,V3)/2, <5 2 - 
a (3, — \/3~) /2, = 0, and a is the C-C bond length. 
We obtain: 
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where c|(fc) and c s (k) are creation and annihilation op- 
erators of conduction (or valance) band with wavevector 
k and spin s, respectively, N is the total number of unit 
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cells in the system, K is a reciprocal unit vector, and 
the summation of ki is over the Brillouin zone. We as- 
sume that inter-band contribution can be ignored, and 
the band index I is dropped for brevity. The interaction 
matrix elements have the form: 

V s3 ' (k^ = ^ [l + (3^1 +^2 -'¥"=3 -^4] S s - S , 
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We focus on the case of lightly doped graphenes with 
kpa <C 1, where kp is the Fermi momentums of the two 
Fermi pockets around the corners of the Brillouin zone 
at ±Qoj as shown in Fig. [TJ We introduce valley index 
A = ± to denote which valley the electron is in, and 
the corresponding valley electron annihilation operator 
is defined as Ck\ s = c s (XQq + k), where k is electron 
momentum relative to the center of the corresponding 
valley. 

It is easy to see that there could exist two forms of 
superconducting pairings in graphene, i.e., inter- valley 
and intra- valley pairings, as shown in Fig. [T] In all the 
previous studies, only the inter- valley pairing, i.e., two 
electrons in a pair have the opposite momentums and re- 
side at the different valleys (Fig.QJa)), is considered. On 
the other hand, the intra- valley pairing, i.e., two elec- 
trons in a pair reside in the same valley, and the two 
valleys in the vicinity of ±Qo act like two independent 
sub-systems (Fig.QJb)), is also possible. To consider both 
possibilities, we enumerate all possible combinations of fcj 
in the proximity of Qo or Qq in Eq. ([3]), and collect all 
the terms relevant to the superconductivity pairings, i.e., 
those terms satisfying either k\ + fc 2 = k% + k± = (for 
the inter- valley pairing) or fci + fc2 = k^ + k4 = ±2Qq (for 
the intra- valley pairing). The effective electron-electron 
interaction can then be written as, 



f7 s < 
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and to the first order of \k\a, \k'\a, the matrix elements 
can be approximated as: 
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where we have made use of approximations: 7*. « 3, 
7±Q +fe = 7±2Q +fe ~ 0, exp(i<p ±Qo+fe ) sa -iexp(=Fi# fc ), 
to the first order of ka, and 9k is the azimuth of k. 



We can then determine all possible pairing states and 
their relative stabilities. First, we consider the case go < 
and g\ > 0. For the inter- valley pairing, we define the 
superconducting gap Ai = -(M/2-/V)X)fc (c-fc-|Cfc+T)> 
and A 2 = -(\g \/2N) £ fc (c-k+i,c k -f) ■ It is easy to 
show that self-consistent gap equations read: Ai = 
A 2 = (|3o|/2)EfeAtanh( ( S J B fe /2)/(2 J B fe ), where E k = 
V( e k - /x) 2 + A 2 , and A = A x + A 2 , e k = ±v\k\ is 
the free electron dispersion of graphene, fi is the chemi- 
cal potential of the system. The equations yield usual 
s-wave pairing with superconducting gap A = Ao = 
2fuvo exp (—l/\go\No) at the zero temperature, where Hujo 
is the cut-off energy of the attractive interaction (e.g., De- 
bye energy in the case of phonon-mediated interaction) 
and No is the density of state per unit cell per valley 
at the Fermi level. We note that the pairing correlation 
function has even parity when exchanging the valley in- 
dexes, and odd parity when exchanging the spin indexes, 
respectively, corresponding to a spin-singlet pairing state 

(| + -) + I - +))(| t4) - I 4-t))- 0n the other hand, Nere- 
is no pairing in the spin triplet channel. 

For the intra-valley pairing, the relevant interaction 
component is U kk , = -(\g \/2) cos(9 k - 9 k >) cxp[iA(6*fc - 
6k')) in Eq. ([5]). We define the superconducting gap 
as A A (fc) = (l/N)J2 k ,U£ k ,(c-kWx t ), and A = ±. 
The self-consistent gap equation reads: A A (fc) = 
-(W E fe < ^ A A (fcQtanh(/3£ Afc 72)/ (2£ Afc Q 

and E xk = \j{tk-Li) 2 + |A A (fc)| 2 . The equa- 
tion has an s-wave solution A A (fc) = Ao and 
a d + id solution A A (fc) = Ao cxp (2iA6>fc), with 
A = 27kJoexp(-4/|< 7o |iVo). 

Second, we consider the case go > and g\ < 0. 
For the inter-valley pairing, the relevant interac- 
tion component is Ukk' = — (3|.9i|/2) cos(0fc — 9k') 
define the superconducting gap 
Ek' u kk' (c- k '- s 'Ck'+ s ). 
gap equation reads, A ss >(k) = 
,(k')t&nh(/3E k ,/2)/(2E k ,) and 
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in Eq. © 

as Ass'(fc) = 
The self-consistent 
-(l/N)Ek' Ukk'A t 

E k = \J (efe - p) 2 + |A ss /(fe)| 2 . The equation 
yields two stable zero temperature solutions in 
p-wave channel: A ss /(fe) = AoCxp(±i#fc) and 
Ao = 2hu>oexp[—A/3\gi\No]. The pairing can be 
either spin singlet or triplet, with the corresponding 
valley states (] + -) - | - +)) and (| + -) + | - +)), 
respectively. [24| 

For the intra-valley pairing, the relevant interac- 
tion component is U kk , = — (3|<7i|/2) exp[iA(0fc — 9k')] 
in Eq. ([5]). We define the superconducting gap 
as A A s ,(fc) = (l/A0£ fe ,^ fe ,<c_ fe , As ,c fe , As >. The 
self-consistent gap equation reads A A S , (fc) = 
-( V^) Ek' Ukk'^Ak') tanh (/3£ fc ,/2)/(2E fc ,) and 

Ek = J (efe — /i) 2 + I A A S , (k) | 2 . The equation yields the 
spin triplet p + ip pairing state at zero-temperature: 
A A S , = A exp(iA6> fc ) and A = 2hu} exp(-2/3\g 1 \N ). 
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The above considerations assume that the Fermi level 
is close enough to the Dirac points, such that the quasi- 
particle dispersion has perfect rotation symmetry within 
each valley. However, real graphene systems only have 
three- fold rotation symmetry. The asymmetry, i.e., 
e±Q +k 7^ e±Q -fcj becomes more pronounced when 
Fermi level is removed from the Dirac points, leading 
to pair breaking that destabilizes the intra-valley pair- 
ing. On the other hand, the inter-valley pairing is not 
affected because the symmetry between tk and e„fc is 
guaranteed by time-reversal symmetry. To see asym- 
metry effect to the intravalley pairings, we expand the 
tight-binding dispersion to the second order, and for the 
valley centering Qq, we have = £Q +k ~ hvk + hk, 
and hk = (hk 2 va/A) sin3#fc. The gap equation for 
the intra-valley pairing when g\ < is modified to 
A^(fc) = -(l/A0E fe < Eft,,A* - f{E k , + h k ,) - 

f{E k ' + h-k')]/(2Ek')- Solving the equation, we find 
that the intravalley pairing states become unstable when 
[i > 2y/ (hv/a)A w 130-\/A/fmeV meV for graphene. 

The discussion can be further extended to the more 
general form of the interaction. We introduce g^p as 
e-e interaction between two electrons resided in the 
unit cell i sublattice a and the unit cell j sublat- 
tice j3 with a, (3 € {A, B}. The Fourier transforma- 
tions are g (q) = (l/N)J2g^ =0 cxp [-iq ■ (Rj - Rj)] and 

ij 

9i(q) = (l/^E^/S CX P H<? ■ (R< - correspond- 

ij 

ing to intra-sublattice and inter-sublattice interactions, 
respectively. In the limit of kpa <C 1, only Fourier com- 
ponents at q = and q = ±2Q are relevant. Moreover, 
we can show that gi(±2Q ) always vanishes because of 
C3 rotational symmetry. The interaction matrix elements 
in Eq. (H} are then modified to: 



Writer fe ) 



ffo(O) 



cos( 



gi(°) r iA(g fc -e,, 
2 

go(0) .9i(0) 

2 2 
ffo(2Qo) 



■ COS(V k 



(8) 
(9) 
(10) 



Following the same procedure, we can extend our anal- 
ysis to the general model. The results are summarized in 
TableHl For the simple model Eq. (JTJ, g (0) = g (2Q ) = 
go, and g\ (0) = 2>g\. Compared to the simple model, the 
general model yields extra possibility of forming inter- 
valley spin triplet s-wave pairing state. 

From Table [H it is straightforward to determine the 
pairing ground states in different parameter regimes. The 
phase diagram is shown in Fig. [2] The most interest- 
ing feature of the phase diagram is the presence of spin- 
triplet intra-valley p + ip pairing state when the inter- 
sublattice attractive interaction g\ (0) dominates, and 
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Table I: The possible superconducting pairing states. The 
momentum dependence of gap function A(fe), the parity in 
exchanging the valley ("valley") and spin ("spin") index, as well 
as the effective pairing strength Vo are shown. The amplitude 
of the zero temperature superconducting gap is related to Vo 
by A = 2hu exp(-l/JV Vo). 



the Fermi level is close to the Dirac points such that 
/i < 2y/ (hv/a)A. It is important to note that the spin- 
triplet intra-valley p + ip pairing state does not break 
the time reversal symmetry (TRS), because the order 
parameters at the two valleys are related by the time re- 
versal, i.e., the order parameter at one valley is p + ip, 
while that at the other valley is p — ip. The origin of the 
intra-valley chiral p + ip pairing can be traced back to 
7r-Berry phase associating with the Dirac point in each 
of the valleys. 



II- 



I3j The topological origin of the un- 
conventional pairing, which is dictated by a band singu- 
larity buried deeply inside the Fermi sea, is different from 
usual mechanisms derived from strong correlations. The 
intra-valley chiral p+ip superconducting state is spatially 
non-uniform, due to the finite center of mass momentum 
(CMM) of Cooper pairs. It is important to note that 
the pairing with the finite CMM does not incur extra 
kinetic energy in graphene, because the electron disper- 
sion does have minimums at ±Qo ( see Fig- HI- This is 
distinct from the another well known non-uniform super- 
conductin g, sta te in conventional systems, i.e., the LOFF 



state, 



14 , Il5j | which incurs large kinetic energy and is 
very fragile. 

The intra-valley triplet p + ip superconducting phase 
in lightly doped graphene, once realized, could be used 
to implement graphene-based topological quantum com- 
putation. This could be achieved if one could gener- 
ate half quantum vortices (HQV) in the superconduct- 
ing phase. [l6| Each vortex core will confine a Majorana 
fermion. The quantum state of a macroscopic system 
with a number of such vortices could be transformed by 
adiabatically moving them, with robustness topologically 
protected by their spatial separations. 



17| 



For graphene, there exists a very useful and unique tool 
for implementing the novel quantum computing scheme, 
i.e., pseudo-magnetic fields generated by strains or topo- 
logical defects. fl8Tj2fl | Using a strain field or topolog- 
ical defect, it is possible to create a spatially localized 
pseudo-magnetic flux, which is not easily achievable using 
a real magnetic field. One can use the pseudo-magnetic 
fluxes as nuclei for creating HQVs. For instance, one can 
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Figure 2: Phase diagram of lightly doped graphene. For given 
go(0) and ,gi(0), the pairing states also depend on the param- 
eters <?o(2Qo) an d Mi alm the corresponding conditions are 
shown for each phase. 



generate a total pseudo-magnetic flux of a half quanta 
within the size of superconducting penetrating length. 
Because of the presence of two independent valleys, it 
will create a pair of HQVs, one for each of the valleys. 
The two Majorana fermions localized inside the vortex 
cores will maintain their independence within time scale 
of inter-valley scattering. [22j On the other hand, to cre- 
ate a single Majorana fermion, one could engineer a total 
pseudo-magnetic flux of a quarter quanta. Because the 
pseudo-magnetic field has the opposite directions for the 
two valleys, a perpendicular real external magnetic field 
will enhance the magnetic field experienced by one val- 
ley while suppressing it in the other. It is then possible 
to generate a HQV only in one of the two valleys. [2f| 
The resulting Majorana fermion will be fully topologi- 
cally protected. Finally, one could carefully design the 
strain field and the matching external real magnetic field 
so that the creation of multiple HQVs is energetically 
favorable, 2l| and the graphene-based topological quan- 
tum computation can then be achieved by tuning the 
strain field. 

To realize the novel superconducting phase predicted 
in this Report, it is important to find a modification 
to graphene systems that introduces an appropriate 
level of carrier density and a sizable effective electron- 
electron attractive interaction from such as electron- 
phonon coupling or electron-plasmon coupling. Q Re- 
cently, it has been predicted by the first principles cal- 
culation that graphene could be transformed supercon- 
ducting by lithium deposition. [Io| While the particular 



system has its electronic structure significantly modified, 
and cannot be expected to support the novel p + ip phase 
we predict, it does give the hope that one may find other 
dopants that serve the purpose. 

In summary, we systematically investigate the super- 
conducting phases for lightly doped graphene. We find 
that an intra-valley spin triplet p + ip superconducting 
phase could be achieved in the certain parameter regime. 
We further show that the novel superconducting phase 
could be used to implement topological quantum com- 
putation in graphene by utilizing pseudo-magnetic fluxes 
created by strain field or topological defects. These pos- 
sibilities make graphene a promising candidate for imple- 
menting topological quantum computation. 
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[23] The other solutions have the form A A (fc) = 
v2Ao cos(6fe + (p) exp (iA^fe), and A w 

1.72fiuo exp (— 4/|go|-^Vo)- The corresponding states 
have the higher energy. 

[24] The degeneracy between the spin singlet and triplet 
states will be lifted in the order of (fc^a) 2 , in favor of 
the singlet state. 



[25] To minimize elastic energy, it is usually necessary to cre- 
ate pseudo- magnetic fluxes in pairs of opposite directions. 
In this case, each pair of the fluxes will create a pair of 
Majorana fermions residing in the different valleys, spa- 
tially separated. 



